home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / ssycon.z / ssycon
Text File  |  1996-03-14  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. SSSSSSSSYYYYCCCCOOOONNNN((((3333FFFF))))                                                          SSSSSSSSYYYYCCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SSYCON - estimate the reciprocal of the condition number (in the 1-norm)
  10.      of a real symmetric matrix A using the factorization A = U*D*U**T or A =
  11.      L*D*L**T computed by SSYTRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE SSYCON( UPLO, N, A, LDA, IPIV, ANORM, RCOND, WORK, IWORK, INFO
  15.                         )
  16.  
  17.          CHARACTER      UPLO
  18.  
  19.          INTEGER        INFO, LDA, N
  20.  
  21.          REAL           ANORM, RCOND
  22.  
  23.          INTEGER        IPIV( * ), IWORK( * )
  24.  
  25.          REAL           A( LDA, * ), WORK( * )
  26.  
  27. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  28.      SSYCON estimates the reciprocal of the condition number (in the 1-norm)
  29.      of a real symmetric matrix A using the factorization A = U*D*U**T or A =
  30.      L*D*L**T computed by SSYTRF.
  31.  
  32.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  33.      condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      UPLO    (input) CHARACTER*1
  38.              Specifies whether the details of the factorization are stored as
  39.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  40.              form is A = U*D*U**T;
  41.              = 'L':  Lower triangular, form is A = L*D*L**T.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) REAL array, dimension (LDA,N)
  47.              The block diagonal matrix D and the multipliers used to obtain
  48.              the factor U or L as computed by SSYTRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      IPIV    (input) INTEGER array, dimension (N)
  54.              Details of the interchanges and the block structure of D as
  55.              determined by SSYTRF.
  56.  
  57.      ANORM   (input) REAL
  58.              The 1-norm of the original matrix A.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSSSSSYYYYCCCCOOOONNNN((((3333FFFF))))                                                          SSSSSSSSYYYYCCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      RCOND   (output) REAL
  75.              The reciprocal of the condition number of the matrix A, computed
  76.              as RCOND = 1/(ANORM * AINVNM), where AINVNM is an estimate of the
  77.              1-norm of inv(A) computed in this routine.
  78.  
  79.      WORK    (workspace) REAL array, dimension (2*N)
  80.  
  81.      IWORK    (workspace) INTEGER array, dimension (N)
  82.  
  83.      INFO    (output) INTEGER
  84.              = 0:  successful exit
  85.              < 0:  if INFO = -i, the i-th argument had an illegal value
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.